SCIENTIFIC 

REPORTS 




SUBJECT AREAS: 

CLIMATE SCIENCES 

SCIENTIFIC DATA 

EARTH AND ENVIRONMENTAL 
SCIENCES 



High temperatures-related elderly 
mortality varied greatly from year to year: 
important information for heat-warning 
systems 



ENVIRONMENTAL SCIENCES Yuming Guo 1 ' 2 , Adrian G. Barnett 1 & Shilu Tong 1 



Received 
12 July 2012 

Accepted 
15 October 2012 

Published 
9 November 2012 



Correspondence and 
requests for materials 
should be addressed to 
Y.G. (guoyuming® 
yahoo.cn) 



'School of Public Health and Social Work, Institute of Health and Biomedical Innovation, Queensland University of Technology, 
Brisbane, OLD 4059, Australia, 2 School of Medicine, The University of Queensland, Brisbane, QLD 4006, Australia. 

We examined the variation in association between high temperatures and elderly mortality (age > 75 years) 
from year to year in 83 US cities between 1987 and 2000. We used a Poisson regression model and 
decomposed the mortality risk for high temperatures into: a "main effect" due to high temperatures using 
lagged non-linear function, and an "added effect" due to consecutive high temperature days. We pooled 
yearly effects across both regional and national levels. The high temperature effects (both main and added 
effects) on elderly mortality varied greatly from year to year. In every city there was at least one year where 
higher temperatures were associated with lower mortality. Years with relatively high heat-related mortality 
were often followed by years with relatively low mortality. These year to year changes have important 
consequences for heat- warning systems and for predictions of heat-related mortality due to climate change. 

There is increasing public health concern for heat-related mortality worldwide as the global climate is 
changing 1 . The Intergovernmental Panel on Climate Change (IPCC) has concluded that without increased 
investments in countermeasures heat waves will cause increased adverse health impacts, including heat- 
related mortality 2 . Urban areas may be particularly vulnerable to heat waves because of high concentrations of 
susceptible population groups and urban heat islands 3,4 . 

Many studies have shown that high temperatures are related to non-accidental deaths 5,6 and cause-specific 
deaths 711 , and other studies have specifically examined the impacts of heat waves 1214 . Gasparrini and Armstrong 
(2011) decomposed the effects of heat into two parts: a "main effect" related to the independent effects of high 
temperatures, and an "added effect" due to heat waves. Results showed that the "main effect" contributed the 
most to the excess risk of mortality, while a smaller added effect appeared during heat waves. Rocklov et al. (2010) 
found an added effect of heat waves with two or more days above the 98 th percentile of temperature. Anderson and 
Bell (2011) examined the modification of heat wave characteristics (e.g., intensity, duration, and timing in 
summer) on mortality. They found higher risks for heat waves that were more intense, longer, or occurred in 
early summer. 

Epidemiological studies of heat-related mortality could be used by decision makers to design an early warning 
system for high temperatures, by establishing a heat threshold and the expected increase in deaths above the 
threshold. Epidemiological studies are also useful for providing a baseline from which to estimate the potential 
health effects of climate change. However, these previous epidemiological studies gave estimates for the whole 
study period, and ignored the potential variability in risks from year to year. Risks may vary from year to year 
because of differences in the at-risk population (e.g., more elderly people), or because of increased adaptation over 
time 15,16 . 

Temperature-related deaths are more pronounced in the elderly as they are more sensitive to temperature 8,17 . 
The thermal regulation system weakens with age; skin sensory perception may diminish and thermal homeostasis 
may decline 18 . This means the elderly are not as physiologically able to get relief from cold and hot temperatures. 
A higher incidence of pre-existing cardiovascular and respiratory disease in the elderly is another reason why they 
are vulnerable to temperature extremes 19 . We therefore restricted our analysis to the elderly as they are a high risk 
group who will have strong associations between high temperatures and death, which gave us the opportunity to 
examine how these associations vary over time and space. 
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Results 

The yearly variability for both the main effect of high temperatures 
and the added effect is shown for each city in Figure 1. In most cities 
there was an extremely wide variation in the yearly estimates, and in 
almost every city both the main and added yearly effects were some- 
times negative. The range in the main effects was greatest in 
Lexington, where the percentage change in yearly elderly mortality 
ranged from —80% to 150%. The smallest range was in San Diego, 
where the percentage change in the added effect ranged from —57% 
to 29%. 

We found a similarly large variation in the associations at the 
regional level (Table 1). In every region both the ranges in the main 
and added effects included zero, meaning that in some years high 
temperatures increased mortality, but in other years high tempera- 
tures decreased mortality. The range in estimates at a national level 
also included zero. 

The estimated percent increase in elderly mortality over time by 
region is in Figure 2. There was a great variation from year to year in 
both the main and added effects of high temperatures. In general, 
years with relatively high increases in mortality were often followed 
by years with relatively low increases, and vice versa. 

According to the Deviance Information Criteria, the Bayesian 
hierarchical model fit was not improved by using a non-linear trend 
over time, so the combined estimates were based on a linear time 
trend. These linear trends by region are in Figure 3. The main effect 
declined over time in all of the regions. The added effects moved 
closer to zero over time in most of the regions. 

Discussion 

This study estimated the yearly risks of elderly mortality associated 
with high temperatures. Both the main and added effects varied 
greatly from year to year in every city and region. In every city and 
region there were years where high temperatures were associated 
with increases in mortality, and other years where high temperatures 
were associated with decreases in mortality. 

Such strong year-to-year variations have important consequences 
for heat-warning systems, which give advanced warning of danger- 
ous temperatures 20 22 . These systems are designed based on historical 
observations between temperature and death in a city, and aim to 
find a temperature threshold above which spikes in deaths occur. 
They are based on average risks over multiple years and do not 
account for the year to year variations in risk shown here. This means 
that these systems are likely to suffer false positives, as in some years 
the anticipated spike in deaths will not occur. They may also suffer 
false negatives, as the spike in deaths could be greater than expected 
or occur at a lower temperature threshold. False positives have fin- 
ancial implications, as much of the preventive actions planned under 
the current system (such as freeing-up hospital beds) are wasted. 
False positives may also undermine confidence in the system if the 
warnings are perceived to be too frequent and without any real 
need 23 . We suggest that people who run heat-health warning systems 
should consider the yearly variability on elderly mortality, and 
update their systems by each summer. 

The great variation in temperature effects shown here also has 
implications for estimating the future health effects of high tempera- 
tures due to climate change 24 . These projections have been based on 
associations between temperature and health that ignore the year to 
year variability in risks. Although the estimated averages from these 
projections are still valid, any confidence intervals will be too narrow 
as they missed the uncertainty demonstrated here. 

There was some visual evidence that years with relatively high 
heat-related mortality were followed by years with relatively low 
mortality (Figure 2). This could be because of changes over time in 
the susceptible pool (or mortality displacement) 15,25,26 . A summer 
with a relatively high level of mortality means less people susceptible 
to heat-related mortality in the following summer. Studies have also 



shown that high winter mortality leads to lower mortality (mortality 
displacement) in the following summer 15,25,26 , as well as high summer 
mortality leading to lower winter mortality 27 . 

The main high temperature effect was generally greater than 
added heat wave effect (Figure 3). The main effect was higher in 
the Northwest and South California, while the added effect was 
higher in Northeast and Upper Midwest. These findings are con- 
sistent with other studies of heat waves and high temperatures 28,29 . 
For example, Anderson and Bell (2011) found that heat waves had 
smaller effects on mortality in the South. 

There was a decline in the main effect of high temperature-related 
death from 1987 to 2000 in all regions (Figure 3). This is consistent 
with previous US studies 7,30 . An increased use of air conditioning 
may be an underlying reason for this reduction in heat- related mor- 
tality 7,30,31 . US studies have shown that associations between high 
temperatures and mortality were related to the level of air condition- 
ing use 28,32,33 , and that the use of air conditioning steadily increased in 
all areas of the United States during 1980 and 2000 7 . 

We defined high temperature days as two or more consecutive 
days with temperatures above the 90 lh percentile during each warm 
season (May 1 to September 30). Using a stricter definition (e.g., 2 or 
more consecutive days with temperatures above the 95 th percentile), 
there were often no high temperature days in some cities. However, 
the results of our study are consistent with a previous study 34 , which 
showed that the main effect of high temperatures is greater than the 
added heat wave effect. 

We did not consider the variability in the effect of temperature 
according to either cumulative effects or by the timing or intensity of 
the high temperatures. Previous studies have shown that individual 
days of extreme heat increased the risk of mortality 35,36 , and that the 
added effect of high temperature on mortality was greater for con- 
secutive high temperature days compared with non-consecutive high 
temperature days 34,37 . Anderson and Bell (2011) found that a heat 
wave's effect on mortality varied by its intensity, duration, and timing 
in the season 37 . 

We did not control for air pollution as a previous study found that 
the association between heat and mortality in the United States was 
independent of air pollution 29 . Other studies have also not controlled 
for air pollution when assessing the impact of temperature on mor- 
tality 25,34,37 . Exposure to air pollution has very different biological 
effects compared with exposure to high temperature, and therefore 
we would expect them to have independent effects. 

The model was fitted well for the data as judged by residuals (see 
Supplementary Fig. SI online). We changed the degrees of freedom 
for smoothing for time, and the degrees of freedom of the splines for 
temperature and lagged temperature. The results were broadly sim- 
ilar (not shown). 

We examined the main and added effects using different high 
temperature definitions (95 th with 2 consecutive days of duration). 
Both the main and added effects of high temperature increased if the 
95 th percentile was used compared with the 90 th percentile (see 
Supplementary Fig. S2 online). 

This study has some strength. The effects of high temperatures on 
elderly mortality are highly variable from year to year throughout the 
US. This finding is important from a public health perspective, 
because it has implications for heat-warning systems and projections 
of the impact of climate change. We used a flexible non-linear model 
to capture the main effects of high temperatures, and consider lagged 
effects up to 10 days. We used a large sample of 83 cities with a wide 
geographic spread across the United States. 

This study has some limitations. We only used 14 years of data. A 
longer period would be able to better examine the changes over time 
in temperature-related mortality. To accurately estimate the auto- 
correlation in the estimates shown in Figure 2 (which we suspect 
would be negative) would require at least 50 years of data 38 . We did 
not have individual exposure data, and so assumed a common 
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Figure 1 | Boxplots of the yearly associations between high temperatures and elderly mortality by city between 1987 and 2000. The main effect is the 
percent change in mortality for an increase in temperature from 75 th percentile to the median of high temperature days. The added effect is the percent 
change in mortality on high temperature days. There are 14 estimates per city. City abbreviations are explained in Supplementary Table SI. 
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Figure 2 | Mean associations between high temperatures and elderly mortality by region between 1987 and 2000 using a univariate meta-analysis. The 

main effect is the percent change in mortality for an increase in temperature from 75 th percentile to the median of high temperature days. The added effect 
is the percent change in mortality on high temperature days. 



exposure to temperature for people in the same city. However, a study 
reported that temperatures from one monitoring site or averaged 
temperature from a network of sites are the same as spatial tempera- 
tures for assessing city-wide temperature effects on mortality 39 . 



In conclusion, both the effects of high temperatures on elderly 
mortality varied greatly from year to year, making it difficult to 
predict the health effects of heat in any one year. There was a gradual 
reduction in high temperature-related elderly mortality. Cities using 



Table 1 | Summary statistics for the estimated yearly high temperature effects on elderly mortality by region between 1 987 and 2000. The 
main effect is the percent change in mortality for an increase in temperature from 75 th percentile to the median of high temperature days. The 
added effect is the percent change in mortality on high temperature days. There are 1 4 estimates for every region. All cells show the percent 
change in mortality 

Main effect Added effect 



Region 
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Industrial Midwest 


-8.1 


0.5 


6.5 


8.9 


16.4 


-6.9 


-5.1 


-2.2 


-0.4 


7.8 


North East 


-16.4 


-1.4 


5.5 


14.2 


19.1 


-9.8 


-5.8 


0.0 


7.0 
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National 
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Figure 3 | Trends over time in the associations between high temperatures and elderly mortality by region using a Bayesian hierarchical model. The 

main effect is the percent change in mortality for an increase in temperature from 75 th percentile to the median of high temperature days. The added effect 
is the percent change in mortality on high temperature days. 



heat warning systems need to be aware of both the long-term trends 
in heat-related mortality and the year to year variations. We recom- 
mend that cities evaluate the performance of their warning systems at 
the end of every summer. Specifically they should check for false 
positives and false negatives, and verify that the temperature thresh- 
old estimated from historical data remains valid. 

Methods 

Data. We used the data from the publicly available National Morbidity and Mortality 
Air Pollution Study (NMMAPS) study 40,41 . This data comprises daily weather, air 
pollution and mortality for 108 cities in the United States from 1987 to 2000. Data on 
maximum and minimum temperatures were from the National Climatic Data Center, 
and daily mortality counts were from the National Center for Health Statistics. Daily 
non-external deaths counts were for residents, excluding injuries and external causes. 
More information is available from the NMMAPS web site (http:// 
www.ihapss.jhsph.edu). 

We excluded small cities (population under 200,000) and cities with more than 
0.5% missing data for air or dew point temperature, which left 83 cities. We stratified 
the cities into seven regions (Industrial Midwest, North East, North West, South East, 
South West, Southern California, and Upper Midwest) (see Supplementary Table SI 
online). We limited analyses to elderly morality (age ^ 75 years) in the warm season 
(1 May-30 September) as we were interested in the effects of heat on a susceptible 
population. Mean daily temperature, calculated as the average of daily maximum and 
minimum temperatures, was used as the main exposure variable. 



Data analysis. A city-specific Poisson regression model was used to examine each 
year's high temperature effect on elderly mortality. In the city-specific model, the heat 
effects were divided into a "main effect" and an "added effect" according to a previous 
study 34 . These two estimates were then combined across cities using a univariate 
meta-analysis to create yearly estimates for each region and for the entire US. The 
city-level estimates were also combined using a Bayesian meta-analysis in order to 
estimate trends over time. 

To examine the city- specific main effect and added heat wave effect in the warm 
season, we used the following Poisson regression model in each city and year: 

Y t ~ Poisson(// ( ), 

/0£(/!,) = a + £r u + uHT,+S(f,3) + ADOW,, (1) 

where Y t is the observed daily death count on day r; a is the intercept; T t j is a matrix 
obtained by applying the distributed lag non-linear model (DLNM) to mean tem- 
perature, p is vector of coefficients for T ti i, and / is the lag days; S(f, 3) is natural cubic 
spline of time, with 3 degrees of freedom for time (1 May-30 September) used to 
control for trends and seasonal patterns in mortality; HT t is a categorical variable for 
high temperature days on day f, HT ( ~0 if day t was a non-high temperature day, 
HT f = 1 if day t was a high temperature day (defined below); DOW f is a categorical 
variable for day of the week, and JL is vector of coefficients. 

We modelled the main effect of high temperatures on day t using the term fiT t j 
which is fitted using a DLNM. A DLNM is a two-dimensional spline which models 
the main effect of heat along dimensions of temperature and lags. DLNMs overcome 
the fact that temperatures within a couple of days are strongly correlated by con- 
straining the lagged (delayed) effects of temperatures using a spline 42 . A cubic 
spline with 5 degrees of freedom was used to model the lagged (delayed) effect of 
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temperature on mortality up to 10 days 34 . A cubic spline with 3 degrees of freedom 
was used to model a non-linear association between temperature and mortality. We 
used the DLNM to give an estimate of the overall temperature effect by summing over 
lag days 43 . This overall estimate was the percent change in mortality between the 
median temperature among high temperature days against the 75th percentile of 
temperature in the warm season. 

The added effect of high temperatures was modelled by «HT t . We defined high 
temperatures as two or more than two consecutive days with temperature above 
the 90 th percentile for each city during each warm season in each city (May 1 to 
September 30). 

To examine the variation in the high temperature effects on mortality from year to 
year, we pooled the yearly national and regional high temperature effects (main effect 
and added effect separately) across cities using a univariate meta-analysis with a 
random effect for each city. We plotted the pooled yearly national and regional main 
effects and added effects from 1987 to 2000. 

We used a Bayesian hierarchical model to estimate trends over time for each region 
and nationally by combin-ing the yearly estimated effects of high temperatures within 
each city and incorporating the estimates' variance 44 . Each city was given a random 
intercept to model its mean heat effect, and a random linear effect of time to model 
linear trends over time. A sensitivity analysis was conducted to examine random non- 
linear effects of time. Model fit was assessed using the Deviance Information 
Criteria 45 . Bayesian hierarchical models have been used to com-bine estimated effects 
across communities in air pollution and temperature studies 7,37 . 

Sensitivity analyses were carried out on the parameters for the city-specific model 
to test the robustness of the results to our assumptions concerning the temperature- 
mortality relationship. We modified the degrees of freedom for smoothing time from 
4 to 8, and the degrees of freedom for the splines for temperature from 4 to 8, and for 
the lag from 3 to 8. 

To test the sensitivity of our results to the definition of high temperatures, an 
alternative high temperature definition was used: 95 th percentile for at least two 
consecutive days. 

The residuals were examined to evaluate the adequacy of the city-specific models. 
The R software (version 2.1 1.0, R Development Core Team 2009) was used to fit city 
specific models. The "dlnm" package of R software was used to create DLNM. The 
"metafor" package of R software was used to fit univariate meta-analyses. WinBUGS 
software (version 1.4) was used to fit the Bayesian hierarchical model. 
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